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Electronic spin and orbital {dd) excitation spectra of (Caa;Lai_a;)(Bai.75_a;Lao.25+a;)Cu30y sam¬ 
ples are measured by resonant inelastic x-ray scattering (RIXS). In this compound, Tc of samples 
with identical hole dopings is strongly affected by the Ca/Ba substitution x due to subtle variations 
in the lattice constants, while crystal symmetry and disorder as measured by line-widths are x in¬ 
dependent. We examine two extreme values of x and two extreme values of hole-doping content y 
corresponding to antiferromagnetic and superconducting states. The x dependence of the spin mode 
energies is approximately the same for both the antiferromagnetic and superconducting samples. 

This clearly demonstrates that RIXS is sensitive to J even in doped samples. A positive correlation 
between the superexchange J and the maximum of Tc at optimal doping is observed. We 

also measured the x dependence of the d^y —>■ d„, 2 _y 2 and d^-^/y^ d„, 2 _y 2 orbital splittings. We 

infer that the effect of the unresolved d^^ 2_^2 d„, 2 ^y 2 excitation on is much smaller than the 

effect of J. There appears to be dispersion in the d^y —>■ d„, 2 _y 2 peak of up to 0.05 eV. Our fitting 
of the peaks furthermore indicates an asymmetric dispersion for the d„, 2 _y 2 excitation. A 

peak at ~0.8 eV is also observed, and attributed to a dd excitation in the chain layer. 

PACS numbers: 74.62.Bf, 74.25.Ha, 75.30.Ds, 78.70.Ck 


I. INTRODUCTION 

Theories built around coupling of the electron spins 
S have become the prominent models for high- Tc 
superconductivity^. A key parameter in these theories 
is the magentic superexchange energy J, which is pre¬ 
dicted to limits or set^i^ the critical temperature for su¬ 
perconductivity. One method of testing this has been 
to compare Tc against J for a variety of cuprates^^— . 
The study of Munoz et al^ resulted in a ATc"*““/A J ~3 
K/meV. However, if the compounds vary in structures 
and nuances, other factors besides J are likely to influ¬ 
ence the J-Tc plot, which are a likely source of scatter in 
the plot of Ref. 5. Another approach has been to mea¬ 
sure the effect of pressure on a single compound. For the 
case of YBCO, Tc has been found to initially increase 
under hydrostatic pressure^r—. Under pressure, J also 
increases^, yielding AT™““/AJ~I.5 K/meV. While sim¬ 
ilar order-of-magnitudes are encouraging, it shows that 
the fluctuations in the slope could be large depending on 
materials or conditions. In fact, Mallet et alX observed a 
negative J-T^ slope in a series of RA 2 Cu 30 y compounds 
with A=(Ba, Sr) i?=(La,..Lu,Y), casting doubt on the 
spin-mediated scenarios. 

Another key parameter thought to strongly affect the 
cuprates is the d^ 22 _y 2 —>• d„: 2 _y 2 orbital splittingH^ii. 
This splitting increases with increasing apical oxygen 
distance dA from the copper-oxygen plane. When the 


splitting grows, it increases the in-plane character of the 
holes, creating a condition favorable for superconductiv¬ 
ity by stabilizing the Zhang-Rice singlet^^ or rounding 
the Fermi surfacei^. A higher dA also reduces screen¬ 
ing from polarizeable charge reservoir layersi^. All three 
options are expected to lead to higher Tc. 

The multitude of different control parameters for T/‘“^ 
emphasizes the importance of measuring their effects in 
isolation. Here we measure both J and orbital splitting 
in (GaxLai_x)(Bai.75_xLao.25-i-x)Gu30y (GLBLCO), us¬ 
ing resonant inelastic x-ray scattering (RIXS). CLBLGO, 
whose phase diagram is shown in Fig. [U is a compound 
which allows the tuning of structural parameters inde¬ 
pendently of the hole doping. Its structure is almost 
identical to YBCO^, but it is tetragonal and its chain 
layers are not ordered. The oxygen content y controls the 
number of doped holes, only slightly affecting the lattice 
parameter. In complimentary fashion, Ga/Ba content x 
changes only structural parameters such as bond length 
a, buckling angles 0, and apical distance dA, while keep¬ 
ing the net valence fixed^i. Additionally, the entire dop¬ 
ing range can be spanned from undoped to overdoped 
for all values of x. Therefore, x tunes both J (through a 
and 9) and orbital splitting (through d^) over the whole 
phase diagram. Moreover, disorder in CLBLGO was 
found to be ^-independent based on the line-widths mea¬ 
sured by techniques ranging from high resolution pow¬ 
der x-ray diffraction^^, Cu, Ca, and O nuclear magnetic 
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resonance^^^— , phonon^l, and ARPES^. 

Intriguingly, both J, as measured in undoped 
CLBLCO samples, and were found to increase 

with X by as much as 40%. In fact, the energy scale of the 
entire phase diagrams, including the magnetic, spin glass, 
and superconducting parts, scale with J— Such scal¬ 
ing was attributed to a superconductivity governed by 
j_27-_30 However, J in the optimally doped samples, and 
the possible effects of the apical oxygen distance, are not 
known. Measuring those are the main objective of this 
work. 

Recent milestones in the technique of resonant inelas¬ 
tic (soft) x-ray scattering (RIXS) have been the mea¬ 
surement of dispersive magnetic excitations in supercon¬ 
ducting cuprates and iron pnictides^^—. This, together 
with single-crystal CLBLCO growths^, have allowed us 
to measure how J varies with x in both underdoped and 
optimally doped CLBLCO samples. The momentum de¬ 
pendence provided by RIXS enables the precise determi¬ 
nation of J based on the spin-wave dispersion. Fortu¬ 
itously, the same probe is also sensitive to the orbital dd 
excitations^!. Here we measure of both effects simulta¬ 
neously on CLBLCO single crystals. We also use x-ray 
absorption spectroscopy (XAS) to verify that the effec¬ 
tive hole dopings are indeed the same when we compare 
families with different x. 

We find that has a positive, but not propor¬ 

tional, correlation with J. The dxy dx 2 -y 2 and 
dxzjyz dx 2 -y 2 splittings also increase with x, but we 
could not precisely isolate the d^^ 2 _x 2 — dx 2 -y 2 excita¬ 
tion. We nevertheless determine that in this system A J 
has a greater contribution (treating it as the indepen¬ 
dent variable) to the change in compared to the 

out-of-plane orbital effect!!^— . We also observed that 
the change in J is very similar in the undoped and doped 
samples, and the slope of the relation is identical 

to that of YBCO under pressure. The RIXS spectra also 
revealed unexpected features, including a peak at 0.8 eV, 
and energy dispersive dd excitations. 

The paper is organized as follows. Section HIl describes 
the experimental details. Presentation and analysis of 
the RIXS data are divided into the sections: pill) RIXS of 
underdoped samples, with focus on the spin excitations; 
m RIXS of doped samples; and m dd, or “crystal 
field” orbital excitations. Section lyTl is a discussion of the 
context of these results, and Section lVlIl is the conclusion. 
Analysis of the doping from the O Al-edge and Cu L-edge 
x-ray absorption spectra is provided in the Appendix. 


II. EXPERIMENTAL DETAILS 

(CaxLai_x)(Bai. 75 _xLao. 25 -i-x)Cu 30 y single crystals 
were grown using the traveling float-zone method!^. For 
each of a;=0.1 and a;=0.4, the under-doped (UD) sam¬ 
ples (in y) were prepared by annealing in argon. The 
near-optimally doped (OD) samples were first annealed 
in flowing oxygen, followed by 100 Atm oxygen pressure 



FIG. 1: (Color Online) Phase diagram of CLBLCO obtained 
from powder samples in Ref. 29, which plots the Neel (Tn), 
spin-glass (Tg), and superconducting (T,) transition tempera¬ 
tures as a function of the stoichiometric oxygen amount y for 
various families {x). The under-doped (UD), medium doped 
(MD), and near-optimally doped (OD) single-crystal samples 
of the present study, are placed in the diagram as triangles, 
black (pointing down) for a; = 0.1 and red (pointing up) for 
X = 0.4. 


for a period of two weeks. The oxygen content for the 
OD samples was confirmed by iodometric titration. The 
oxygen content for UD sample was set based on the pro¬ 
cedures used for powders so as to be in the antiferromag- 
netically long-range ordered phase!^. The T^s of the 
X = 0.4 and a; = 0.1 OD samples were measured by mag¬ 
netic susceptibility using a SQUID magnometer. They 
were 78 K, and 46 K respectively, indicating a close to 
optimal doping condition. The place of the samples in 
the phase diagram is shown in Fig. [TJ 

Soft x-ray absorption spectroscopy (XAS) and 
RIXS measurements were conducted at the ADRESS 
beamline!! at the Swiss Light Source of the Paul Scherrer 
Institut. The sample environment was ~10 K in vacuum. 
The sample surfaces were cleaved c-axis faces, mounted 
such that the a- (or equivalent b-) axis was in the horizon¬ 
tal scattering plane. For XAS, to obtain incident polar¬ 
ization approximately parallel to the c-axis, the sample 
was rotated to 10° from the grazing incidence condition. 
Refer to the Appendix for the detailed XAS results. 

RIXS spectra were measured in the horizontal scatter¬ 
ing plane. Measurements were done for both horizontally 
and vertically polarized incident beams, corresponding to 
TT and a polarizations respectively. The incident energy 
was set to the first main peak in the Cu Lui XAS at 
932 eV. The detector was fixed such that the two-theta 
scattering angle with respect to the incident beam was 
130°. Throughout this article, we refer to the in-plane 
momentum transfer q in reciprocal lattice units of 2t: j a, 
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where a is the lattice constant of the crystal. We de¬ 
fine q as the component of the total momentum change 
of the photon which is parallel to the sample ab planed. 
Our sign convention is grazing incidence corresponding to 
negative q. The variation of a with each is accounted 
for in calculating q, but is not significant on the g-scale. 
In our scattering configuration q is always along the (1 
0 0) direction, and its magnitude is changed by rotat¬ 
ing the sample away from specular reflection. Therefore 
the total momentum transfer is Q—{q, 0, L) in tetrago¬ 
nal notation. The grazing incidence condition was used 
to calibrate the q position. This calibration was found 
to be valid by measuring E ys. q dispersion for both 
positive and negative q (see Section nni). 


III. MAGNONS IN UNDERDOPED SAMPLES 

Typical RIXS spectra for the UD samples are com¬ 
pared in Fig. for a; = 0.1 and x = 0.4. The vari¬ 
ous panels of Fig. [2] zoom in on different energy scales. 
The intensities are normalized to match at the strong dd 
peak in Fig. HKb), and the energies are shifted so that 
the quasielastic peaks in Fig. [UJc) are centered at zero. 
Fig. UKa) shows the relatively high-energy part of the 
spectra. First, we note that the x = 0.1 and x = 0.4 tails 
going down from 5 eV overlap closely. Secondly, there 
is a peak at around 4.5 eV which is in the energy range 
of charge transfer excitations across the Hubbard gap. 
The feature is shifted to higher energy for x = 0.4, and is 
also present in the doped samples. A doping-independent 
feature at similar energy was studied in LSCO with Cu 
AT-edge RIXS^. Fig. HKb) is an overall view of the spec¬ 
tra including both the intense peak encompassing the dd 
excitations between 1.5 and 1.8 eV, and the lower energy 
peaks, which are much lower intensity, but still visible on 
this scale. Comparison of the x = 0.4 and a; = 0.1 spec¬ 
tra over a broad range reveals that generally, most of the 
excitations are at slightly higher energy for x = 0.4. This 
increased energy is ubiquitous both for the magnon exci¬ 
tations covered in this and the following section, and for 
the dd excitations in Section El We note that both the 
high energy tails in Fig. [IJa), and the quasielastic peaks 
in Fig. He), are aligned in energy for the two samples. 
We will later show that the magnon and dd excitation 
energy increases can be directly attributed to the change 
in lattice parameters. 

Fig. He) zooms in on the low energy range. At zero 
energy is a quasielastic peak, which depends on a com¬ 
bination of finite g-resolution of the instrument and the 
crystal mosaic of the sample. In our analysis, the energy 
scale of each spectra is shifted according to the center en¬ 
ergy of the quasielastic peaks. In similar measurements 
done by Braicovich et al^ for La 2 Cu 04 , in which the 
quasielastic peak was much lower than it is here, a fea¬ 
ture at around 80-90 meV was observed, with about a 
fifth of the magnon intensity. This was attributed to a 
resonantly enhanced optical phonon. We did not detect 



FIG. 2: (Color Online) The main features of the typical 
CLBLCO RIXS spectra for x=0.1 UD (black) and x=0A UD 
(red) at various energy ranges; (a) 3-8 eV (b) 0-3 eV and (c) 
below 1 eV. In this example, the incident beam is rr-polarized, 
and (jt=-0.34 


such a phonon and it is not included it in our analysis. 

The peak associated with magnons is found in the 0.2- 
0.4 eV range of Fig. He). Comparison of the data for 
x = 0.4 (red) and a; = 0.1 (black) clearly shows that the 
X = 0.4 peak is shifted to higher energy. Thus the main 
result that J is higher in the a; = 0.4 sample than in 
the a; = 0.1 samples is clearly evident already in the raw 
data. 

There is also a peak at 0.8 eV. Its intensity is highest 
(comparable to the magnetic peak) at negative q for tt po¬ 
larized scattering, but can be seen elsewhere (see Fig. H 
and is always stronger for the x = 0.4 sample. Where 
it is large, it was incorporated into our fitting for the 
magnons, described below. It has only slight dispersion 
of < 0.05 eV, unlike the new mode recently observed by 
Lee et alA^. It would be surprising if the ^ 0.8 eV peaks 
were one of the three dd excitations, which are expected 
to be above 1.5 eV— . On the other hand, a dd excitation 
in the chain layer would be more plausible. Since half 
of the non-apical oxygen ligands around each Cu atom 
are missing in the chain, the Coulomb energy cost for a 
chain dd excitation should also be about half of a plane 
dd excitation, which corresponds to this ^ 0.8 eV. 

A sample pair of spectra corresponding to the a and tt 
polarizations at the same q are shown in Fig. Ha) and (b). 
To extract the magnon energies, fitting was done over 
the range shown in Fig. H Each spectrum was modeled 
as a sum of quasielastic peak, magnon, 2-magnon, with 
(when visible) an additional peak at 0.8 eV, and a tail 
from the dd excitations. The spectral weight of the 2- 
magnon relative to the magnon is generally different for 
the a and tt polarizations, resulting in a shift in the peak 
energy for the different polarizations. As in Ref. 33, the 
fitting is done for both polarizations simultaneously. The 
energies and widths of the magnon and 2-magnon peaks 
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FIG. 3: (Color Online) The RIXS spectra in underdoped sam¬ 
ples at 5=0.375 for (a) a and (b) tt polarized configurations, 
corresponding to vertical and horizontal polarizations of the 
incident beam for scattering in the horizontal plane. The 
spectra of x=0.1 (black) and a;=0.4 (red) samples are com¬ 
pared. The magnon (~0.25-0.3 eV) and 2-magnon (-^0.4 eV) 
component peaks are indicated as solid and dashed lines, ob¬ 
tained from simultaneous fitting of both (a) and (b) spectra. 
The dotted line shows a fit of the tail of the higher energy 
dd excitation. A peak around at 0.8 eV, much stronger in 
the a:=0.4 sample, is also present. The total fits are shown as 
solid lines crossing the data. 


were constrained to be the same for both polarizations, 
as indicated by the vertical lines in Fig. [31 The lineshapes 
as a function of energy v used for all of the excitations 
was a damped harmonic oscillator response in the form of 
a Lorentzian, weighted according to “detailed balance” : 

^ ( (r/2)^ _ (r/2)^ \ 

\{y-VRY + {Vl2Y (i. +1.^)2+ (r/2)2; 

( 1 ) 

where T is the sample temperature, /cs is Boltzmann’s 
constant and fit parameters vr and F are the energy 
and intrinsic width respectively. Each S{v) was then 
convolved with a Gaussian representing the resolution 
function of the spectrometer, to produce the components 


shown in Fig. [S] The fits for all spectra (more than 
80) were excellent and are shown in the Supplementary 
materials^. The dispersion of vr for the magnon compo¬ 
nents is plotted in Fig 01 The horizontal q-axis for each 
sample was corrected by a slight shift (0.013 for x=0.1 
and 0.022 for x=0.4) to make each dispersion symmetri¬ 
cal about the origin. 

The dispersions are fit to a theoretical expression for 
acoustic-mode dispersion in the double layer cuprate 
YBCOi^: 


E = 2 J • ( 1 - 7 ^( 5 ) + ( Jx/2 J) • (1 - 7 ( 5 )) )" (2) 

for in-plane magnetic exchange J, with interplane cou¬ 
pling J± set to 15 meV, and 7 ( 5 ) = 0.5 • (cos(27r ■q) + l). 
There is also in principle an optical mode^, but it resides 
quite close to the acoustic mode, except at low q, where 
the errorbars are high. In our fitting, J± was fixed at 15 
meV (which is similar to YBCO^), so the only free pa¬ 
rameter was J. The fits are shown as the lines in Fig. 01 
Eq. 01 captures the non-linearity of the data, particularly 
well on the negative q side. The resultant J values were 
134 ± 1 meV for x = 0.4 and 120 ± 1 for x = 0.1. 

These values should be compared with detailed ab 
initio calculations done by Petit and Lepetit for opti¬ 
mally doped CLBCO^. Those yielded mean values of 
J = 132 meV for x = 0.4 and J = 110 meV for x = 0.1. 
The X = 0.4 results are in excellent agreement between 
theory and experiment, while there is a 10 meV differ¬ 
ence for X = 0.1. We show in the next section that the 
dispersion of the UD and OD samples are similar. 



FIG. 4: (Color Online) The dispersion along the (1 0 0) di¬ 
rection, of the fitted energies ur of the single-magnon compo¬ 
nents of the x=0.1 (black) and x=0.4 (red) UD samples. Fits 
to theoretical acoustic magnon dispersions of Ref. 43 for free 
parameter J, with J± fixed to 15 meV, are shown as black 
(x = 0.1) and red (x = 0.4) lines. 

The other fit parameters are plotted in Fig. |SDa)-(d). 
In Fig. jSKa), the 2-magnon energies at low q are close to 
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0.3 eV. This magnitude is within the range of the recent 
2-magnon Raman study in this material by Wulferding et 
al^, who measured energies of 0.29 — 0.35 eV in various 
samples. In addition, the sign and magnitude of disper¬ 
sion of the 2-magnon of about 0.1 eV in Fig. [5ja) is rea¬ 
sonably consistent with the ^80 meV measured with O 
K-edge RIXS by Bisogni et alA^ in La 2 Cu 04 (see Figure 
6 of Ref. 44). Fig. [SKb) plots the ratio of the intensities 
of the 2-magnon to the 1-magnon components. They fall 
on the same curve for x = 0.1 and x = 0.4, which is 
expected since the excitations in both should have the 
same symmetries. 

Fig. [a c) and (d) show the intrinsic widths F for the 
magnon and 2-magnon, which are 100-150 meV and 
> 300 meV respectively. These are wider than expected. 
The 2-magnon widths observed in the Raman study^l 
were only ^100 meV, while the magnon width is expected 
to be resolution limited on this scale. It is not clear if 
the large width originates in the fitting or sample. There 
is some intrinsic disorder in the site occupation between 
Ca, Ba, and La atoms, which could be a potential cause 
of an intrinsic magnon width. But if so, we note that 
the widths of a;=0.1 and a;=0.4 are about the same, in¬ 
dicating that X does not affect disorder. Nevertheless, 
our analysis: (1) fit all of the data excellently with min¬ 
imal number of parameters, (2) resulted in a realistic 
dispersion curve with J values which are in good agree¬ 
ment with Ref. 45, and (3) 2-magnon energies at low q 
are consistent with the 2-magnon energies measured with 
Raman scattering^!, and (4) 2-magnon dispersion is con¬ 
sistent with O RT-edge value for La 2 Cu 04 from Ref. 44. 



FIG. 5: (Color Online) The other fit parameters for the x=0.1 
UD (black) and a;=0.4 UD (red) samples: (a) 2-magnon en¬ 
ergy (b) ratio of 2-magnon to magnon intensities in the tt- 
polarized spectra (c) intrinsic magnon width (FWHM) and 
(d) intrinsic 2-magnon width. 


IV. PARAMAGNONS OF OPTIMALLY DOPED 
SAMPLES 

Here we estimate the change in J in the supercon¬ 
ducting samples. For doped cuprates, Le Tacon et al. 
found that the lifetime broadening of the spin excita¬ 
tions make the widths too broad to distinguish between 
magnon and 2-magnon, and instead they are replaced 
by a single “paramagnon” peak!!. A typical spectrum 
for the OD CLBLCO samples is shown in Fig. |6l As 
in Ref. 33 we replaced the magnon and 2-magnon with 
a single magnetic component, retaining the lineshape of 
Eq. [T] Only the elastic intensity, “paramagnon” peak, 
and tail from the dd were included in the fits. Most of 
the q’s measured were positive, and there were no strong 
0.8 eV peaks. Since the peak position is generally dif¬ 
ferent for TT and a polarizations, due to different weights 
of the 2-magnon and magnon contributions (as seen in 
Fig. [3]), both could not be fit simultaneously with one 
peak. We therefore chose to use only the tt polarization. 
The single peak of Eq. [T] plus background fit quite well 
to the data; fits to all of the spectra are shown in the 
Supplementary materials^. As seen in both the fits and 
raw data of Fig. [51 the paramagnon for x = 0.4 is shifted 
with respect to x = 0.1 and extends to higher energy, 
which was also typical for the other q’s. We note that 
in Fig. |6] the dd tails from high energy are the same for 
X = 0.1 and x = 0.4. 



FIG. 6: (Golor Online) 7r-polarized spectra of a;=0.1 (black) 
and a;=0.4 (red) OD samples, at q=0.266. The fit of the 
magnetic component to an asymmetric Lorentzian plus back¬ 
ground, and the magnetic peak component itself, are shown 
as lines with the same color code. The high-energy dd tail of 
a;=0.4 closely coincides with that of a;=0.1, and is shown as a 
dotted line. 

A series of spectra for progressively higher q are plotted 
in Fig. [3 for (a) the UD samples and (b) the OD sam¬ 
ples. The spectra in Fig. [71(a) for the UD samples were 
obtained by subtracting all of the fitted components (see 
Section IITTI) from the raw data, save for the magnon and 
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2-magnon contributions. The same procedure is applied 
to the OD samples in Fig. [TJ^b), by subtracting the non¬ 
magnetic contribution. The q positions are similar for 
Fig.EJa) and Fig.[71[b). Both pairs of spectra in Fig.[71Ja) 
and (b) are centered below 0.2 eV at low q (bottom spec¬ 
tra), and by g=0.4 (top spectra) they dispersed to 0.3 
eV. This similarity suggests that the J comparison for 
the UD spectra, which is generally easier to precisely de¬ 
termine, is also valid for the superconducting case. It 
also would seem to argue against the scenario of intra¬ 
band excitations (as opposed to paramagnons) which was 
recently proposed by Benjamin et al4^, since the OD and 
UD spectra have the same energies. 



Energy Transfer (eV) 



Energy Transfer (eV) 


FIG. 7: (Color Online) Comparison of the background- 
subtracted spectral intensity {Iba) between (a) UD and (b) 
OD samples measured at similar q positions. The spectra 
were obtained by subtracting the quasielastic, dd, and 0.8 eV 
fitted components (if present) from the 7r-polarized spectra. 
The dashed vertical lines are guides for the eye. 


The value of J cannot directly be determined from 
the “paramagnon” spectra. The fitted energy parame¬ 
ter un of the asymmetric lineshape in Eq. [T] does not 
have the same well-dehned meaning as in the two-peak, 
two-polarization fits used in section EH This is be¬ 
cause the peak htted-for here encompasses both magnon 
and 2-magnon components, weighted by some unknown 
amount depending on scattering cross-section for each 
(one can refer to Fig. EKb) for the UD case). In¬ 
stead, for comparison purposes we use the center-of-mass, 
namely, the statistical mean energy {Em ) = ^ 
of the background-subtracted magnetic spectra lbs (E) of 
Fig. Eb). While this definition is arbitrary, for a given 
q, it should be roughly proportional to J for any two 
samples, since both magnon and 2-magnon energies are 
proportional to J. 

(Em) is plotted as a function of q in Fig. |8] for x = O.I 
(black circles) and x = 0.4 (red squares). For all but 
the last, i t is h igher for x = 0.4. The average over these 
q points, {Em), is 0.33 eV for x = O.I and 0.36 eV for 


X = 0.4. Assuming proportionality, we interpret this as 
a 9% increase in J from x = Q.\ to x = 0.4. By compar¬ 
ison, the percentage increase for the (more precisely de¬ 
termined) J’s of the UD samples in Section ED is 11.7%. 
Considering the broad widths of the OD spectra and the 
somewhat cruder method of estimating their AJ, this 
estimated increase is quite close to the UD case. 

In the inset of Fig. [5] we present the negative magneti¬ 
zation measurements of the two superconducting samples 
used for RIXS. There is a clear difference in their T^. The 
main observation of this work is that the sample with 
higher Tc also has higher J. It was also demonstrated 
here that RIXS can distinguish samples with small dif¬ 
ferences in J even in the optimally doped case. 



FIG. 8: (Color Online) The dispersion of the energy center- 
of-mass {Em) of the magnetic peak of the OD samples as 
described in the text. To emphasize their Tc variations, the 
inset shows the magnetization versus temperature of the two 
samples, normalized to their maximum diamagnetic responses 
at low temperature. 


V. CRYSTAL FIELD {dd) EXCITATIONS 

The dd spectra of our UD samples were generally 
sharper than for our OD samples, so we focus on the 
former. The dd excitation spectra of the UD samples 
are plotted in Fig. IHlfor selected q’s for the a:=0.1 sam¬ 
ple. The spectra of the a:=0.4 sample was qualitatively 
similar in the main features, but with slightly higher en¬ 
ergies (see Fig. H^b)). All of the spectra and fittings 
for the full range of g’s are presented in the Supplemen¬ 
tary materials^. The centering of the quasi-elastic peaks 
of all of the spectra are also shown in the Supplemen¬ 
tary materials to be accurate within '^lO meV. At least 
two peaks are clearly resolved, at eV and ^1.7 eV, 
with the intensity of the 1.7 eV peak becoming relatively 
stronger with increased q. We fit the tt and a polarized 
spectra simultaneously to a sum of Gaussians, constrain¬ 
ing the parameters of widths and energies to be the same 
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FIG. 9: (Color Online) dd spectra of the UD sample for a;=0.1 
at representative q positions for (a) a-polarization and (b) 
TT polarization. The three Gaussian components are indi¬ 
cated as dashed/dotted lines, and the total fit as solid lines. 
The filled black squares/dashed lines correspond to positive q, 
the empty blue circles/dotted lines correspond to negative q 
(closer to grazing incidence). The vertical lines and red circles 
are guides for the eye. 


for both polarizations. Three Gaussians worked best. 
They are shown in Fig. 121 The zero energies are defined 
by the elastic peaks (see Section Hill) . As can be seen in 
Fig. 12] the widths successively increased from the low to 
high energy peaks. 

To assign the peaks, we refer to the work of Sala et 
alM, who studied dd excitations with Cu L-edge RIXS 
in a variety of cuprates. They found excellent agreement 
between the observed polarization and q dependence, and 
their cross-section calculations. The compound studied 
in that work which is structurally similar to CLBLCO 
is the double-layer 123-cuprate NdBa 2 Cu 307 (NBCO). 
In what follows, E^y, and E^^ 2_^2 refer to 

the energies of the orbital transitions dxy —f d^^-y^, 
dxz/yz dx 2 -y 2 , and —>■ dx^-y^ respectively. 

The NBCO spectra had two prominent peaks at 1.52 eV 
and 1.75 eV, which the authors of Ref. 42 assigned to 
Exy and Exz/yz- E^z^_x 2 was calculated to be 1.97 eV, 
but it was not visible in their spectra. As q increased, the 
cross-section of the 1.75 eV peak increased relative to the 
1.5 eV peak. These results, both the energies and cross- 
section g-dependence are very close to what we observe 
for CLBLCO in Fig. |21 We therefore likewise assign the 
1.5 eV peak to Exy and the 1.7 eV peak to Exz/yz- Fur¬ 
thermore, the energy of the broad third Gaussian com¬ 
ponent in Fig. 12 ] happened to lie very close to 2 eV, with 
zone-averages (standard deviations) of 1.97(0.03) eV, and 
2.00(0.1) eV, for a;=0.1 and x=QA respectively. While 
this energy is in excellent agreement with calculations 
for E^z'^-r'^ in NBCO^, and for YBCO^ the broadness 
makes it difficult to identify with certainty. 



q { 2 n/a) 


FIG. 10: (Color Online) The fitted dd energies Exy and Exz/yz 
which resulted from the simultaneous fitting to both tt and 
a polarizations, plotted for each q position. The energies of 
a;=0.1 (black) and 2;=0.4 (red) are plotted. 


The g-dependence of Exy and Exz/yz are plotted in 
Fig.|Tn]for x=Q.l and x=0.4. Surprisingly, there appears 
to be some dispersion in the energies. The dxy excitation 
for a:=0.1 shows flat dispersion near the zone center, up to 
around |q|=0.2, but beyond this it exhibits negative dis¬ 
persion of the order of 0.05 eV towards the zone bound¬ 
ary. This Exy dispersion is quite symmetrical about g=0, 
up to g=0.35. The dispersion could also be seen from the 
raw data. In Fig. | 2 ] the red circles mark the low-energy 
peaks of the negative q branch. They also mark the posi¬ 
tive branch, but it is harder to see for high |q|, especially 
for cr-polarization. For the bottom two spectra, corre¬ 
sponding to |(7|=0.09 and \q\=Q.22, the peaks are aligned 
with the vertical dashed line, but by |q|=0.38, the peak 
of the raw data is visibly shifted to the right of the line 
by about 50 meV. A non-zero dispersion suggests propa¬ 
gation of the orbital excitation. For x = 0.4, the dxy ex¬ 
citation shows similar dispersion on the negative branch, 
but its magnitude is roughly halved. Exz/yz also shows 
dispersion, but is not symmetrical about <7=0. The fitted 
energies for a:=0.4 especially show a linear trend with a 
dispersion of almost 0.1 eV between q = ±0.2. Unlike 













































the Exy dispersion, the E^^z/yz dispersion is not obvious 
from the raw data itself due to the wider peaks, and only 
becomes apparent after the fittings. Although it seems 
counter-intuitive, asymmetric dispersion may happen in 
the presence of spin-orbit interaction. It has already been 
observed in the spin-wave of Fe ultrathin films^, for ex¬ 
ample. But, as far as we know this would be the first 
observation of asymmetry in the dispersion of a dd exci¬ 
tation. 

We can check whether E^y scales properly with the 
lattice parameter. As pointed out by Sala et E^y oc 
a“”. Averaging the energies of Fig. (TU] over the zone 
yields, for x = 0.1 {x — 0.4), E^y =1.46 (1.52) eV and 
Exz/yz =1.69 (1.75) eV. The corresponding a values for 
X = 0.1 is a=3.91 A and for x = 0.4 is a=3.88 A. This 
yields n=5.1 remarkably close to the theoretical single¬ 
ion crystal field model’s value of n=5. 


VI. DISCUSSION 

Analysis of the UD spectra in Section Hill provided ex¬ 
plicit J values of 120 meV (x=0.1) and 134 meV (a:=0.4). 
The corresponding for these x values are 57 K and 
80 K respectively^^. In section IIVI we found that the 
change in J for doped samples is comparable to the un¬ 
doped case, and the two dopings furthermore exhibit very 
similar dispersions of the spin-excitation spectra (refer to 
Fig. El). It is therefore justified to apply the UD values of 
J to the superconducting case, as has been assumed to 
be valid in other works^j^l. With x as an implicit param¬ 
eter we find that AT™“’'/AJ=1.64 K/meV. This is the 
same order of magnitude of the average slope obtained 
from the study of Munoz et al<^ of several cuprates hav¬ 
ing different numbers of layers (~ 3.2 K/meV) . It is even 
more closely aligned with the initial slope for YBCO un¬ 
der hydrostatic pressure (^ 1.5 K/meV}i. Moreover, the 
increase of J of 11.7% from x=0.1 to x=QA determined 
for the UD samples in Section m is in close agreement 
with the estimation of 11.9% we obtain by using a simple 
J a cos^ 9/a^'^ rule^S. In addition, E^y scales as expected 
with distances. These results indicate that the in-plane 
energies J(x) and Exy(x) depend purely on in-plane pa¬ 
rameters, without secondary effects arising from different 
Ca/Ba ratios. We speculate that the d^z^-r"^ d^^-y^ 

peak, which we could not properly resolve, behaves as 
expected from the lattice parameters variations between 
different CLBLCO families. 

Whether T^°‘^{x) likewise depends only on the in¬ 
plane parameters is not a priori clear, since the out-of- 
plane lattice parameter c, and apical oxygen distance dA 
are also functions of x. In fact, a number of studiesi^^— 
focused on the effect of dA and E^^ 2 _j .2 on in var¬ 

ious cuprate systems. We now assess the relative impor¬ 
tance that these have for T^°‘^{x). 

Since YBCO and CLBCO share very similar struc¬ 
ture and lattice parameters, it is relevant to compare 
the two. The values of AT^°'^/AJ observed in the 


pressure-dependence of YBCO on one hand, and in the 
x-dependence of CLBLCO observed here on the other, 
are very similar. Hydrostatic pressure compresses the 
c-axis, decreasing the apical oxygen distance dA and in¬ 
creasing Tc- In contrast, when increasing x (and T™““) in 
CLBLCO, dA increases^. That AT^°‘^/AJ is the same 
for YBCO and CLBLCO, in spite of dA changing in the 
opposite sense, leads us to conclude that dA variations 
do not play a major role here in determining . 

Another way to reach this conclusion for CLBLCO is 
to estimate the effect of the change in dA on by 

comparing with other studies. A sensitivity of roughly 
^30 K/A, was shown across various cuprates by 
Johnston et ah— (see Figure 1 of Ref. 49). In CLBLCO 
powder^ as x increases from 0.1 to 0.4, dA increases by 
^0.05 A^S. On that basis, the effect of AdA on Tc in 
CLBLCO would be less than 2 K. 

A similar effect of AdA on Tc results from the theoret¬ 
ical calculations of E^z'^-r^idA) by Sakakibara et alJ^. 
They calculated the Eliashberg eigenvalue A which sets 
a limit on T^. From their calculations, an upper limit of 
<125 K/A can be set, which is still too small to 
account for the Tc variations in CLBLCO. 

Taken together, the above comparisons suggest that 
AdA in CLBLCO has very little impact on By 

eliminating this out-of-plane influence, it becomes more 
likely that the change in observed between differ¬ 

ent families of CLBLCO is due to variations in J. While 
rpmax increases by 40% (Fig. 1) from x=0.1 to x=0.4, J 
as measured by RIXS only increases by ~ 11.7%. From 
other methods, the corresponding increase in J for sam¬ 
ples with the same in-plane hole underdoping was deter¬ 
mined to be: 21% from the 2-magnon Raman peaks^l, 
26% from angle resolved photoemission spectroscopy^, 
20% in ab initio calculations^, and 40% by /xSR with 
extraction of J from Tf^. With the exception of the 
latter, these estimates were all considerably less than the 
increase in This suggests that the J dependence 

of Tc is not proportional, as predicted by some exchange- 
driven theories^. If a linear relationship extends down 
to it would imply a threshold J for supercon¬ 

ductivity. 

VII. CONCLUSION 

To review, we measured the O AT-edge and Cu L-edge 
XAS, and RIXS spectra at the Cu L-edge, in both under¬ 
doped and optimally doped CLBLCO single crystals of 
X = 0.1 and x = 0.4 families which have different 

From the electronic structure of the XAS spectra, sim¬ 
ilar hole dopings in the superconducting samples of the 
different families were confirmed. As it turns out, doping 
does not have a critical effect on the magnon dispersion, 
besides a broadening of the peaks. The relative change 
in magnetic energies between a;=0.1 and a;=0.4 are fur¬ 
thermore similar for the doped and undoped cases. This 
demonstrates that RIXS can distinguish between samples 
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of slightly different J even in the doped case. 

The main dd excitations were also examined and un¬ 
expectedly dispersion of up to 0.05 eV was observed, 
raising the possibility that these orbital excitations can 
propagate. More intriguingly, the dispersion of the exci¬ 
tation from the d^z/yz orbit appeared to be asymmetric 
about (7=0. Higher resolution studies would be needed 
to clarify this dispersion. In the UD samples, an addi¬ 
tional 0.8 eV peak was observed, and attributed to a dd 
excitation in the chain layer. 

Finally, there is a positive correlation between 
and J with a slope consistent with the pressure depen¬ 
dence of both parameters in YBCO. The measured spin- 
wave energies change with x by an amount that would be 
expected from purely in-plane lattice constants change. 
Furthermore, it is concluded that the apical oxygen dis¬ 
tance does not change enough with x to have a signifi¬ 
cant effect on These points suggest that the 

variation with x in CLBLCO is purely an in-plane effect 
driven by orbital overlaps. 
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Appendix: XAS analysis 

In addition to determining the resonance energy 
needed for RIXS, XAS also provides valuable informa¬ 
tion about the number of holes present in our samples. 
We measured the XAS of the single crystal a:=0.4 OD 
and UD, and a:=0.1 OD, and UD samples. In addition, 
an a;=0.1 sample of intermediate doping (MD), estimated 
to be just before the onset of superconductivity, was mea¬ 
sured. We used much of the same approach for analysis as 
was used by Agrestini et for treatment of GLBLCO 
powder. The clearest and most systematic spectra were 
at the Gu L/jj-edge when the electric field is polarized 
along the c-axis (with a 10° misalignment from axis), 
and at the O AT-edge absorption when the electric field 
is polarized parallel to the a6-plane. These are shown 


in Fig. [TT] and Fig. [12] respectively, after subtracting a 
background in the form of an inverse tangent function as 
shown in the insets. Referring to Fig. [TTJ the data were 
normalized so as to have the same maxima of peak A for 
all samples, which comes from the Cu 3d® —>■ Gu 2p3d^® 
transitio n®®!^® . 

The low energy edges of the A peaks of all samples 
match perfectly, with the exception of a:=0.4 UD whose A 
peak is shifted to slightly lower energy. The second peak 
B is at the same energy for all samples. It corresponds 
to the same absorption process as A, but in the presence 
of a ligand hole, namely Cu 3d®Z —?> Cu 2p3d^®jj ®®i^® . 
It is clear that peak B becomes less and less intense as 
the doping decreases, but is roughly the same between 
a:=0.1 and x=0.4 for identical nominal dopings. A third 
peak C appears for the UD samples ~3 eV from peak A. 
It is quite strong for x=0.1 but is only a small bump for 
x=0.4. Such a peak is associated with charge transfer 
excitations to the upper Hubbard band^i^. A satellite 
peak around that energy has been related to the chain 
layer in the 123-compounds^. 



FIG. 11: (Color Online) The x-ray absorption spectra at the 
copper Liii-edge, after background subtraction, for the four 
CLBLCO samples; x=0.4 and 2;=0.1 optimally (OD) doped 
and underdoped (UD) samples, and an additional a;=0.1 sam¬ 
ple at medium doping (MD). The electric field was aligned 10° 
from the c-axis. The three main peaks are labeled A and B 
and C. The inset shows an example fitting of the background. 

The number of holes can be determined from the rel¬ 
ative B peak intensit y®®'^^ . The spectra were fitted 
to three Lorentzians, as shown in Fig. [T3| for the OD 
samples. The ratio of the areas of the components, 
H/(A -I- H), for OD x = 0.4 and x = 0.1 samples were 
0.652±0.01 and 0.657±0.01 respectively, indicating iden¬ 
tical hole doping for both samples. Additionally, we can 
estimate y and the total number of holes in a unit cell 
including chains and planes, h. Roughly 20% of h is ex¬ 
pected to be in each planed. From the measured of 
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FIG. 12: (Color Online) XAS at the oxygen K-edge after 
background subtraction, with electric field parallel to the ab 
plane. The data were normalized to have the same maximum 
intensity. 


TABLE I: Table of parameters determined from the Cu L- 
edge absorption spectra. The columns are sample, relative 
area of the B peak, estimated total number of holes h, and 
estimated oxygen content y. As described in the text, for 
the first two rows y was estimated based on Tc, and then h 
calculated. In subsequent rows h was calculated first, followed 
by V- 


Sample 

B/{A + B) 

h 

y 

a:=0.4 OD 

0.652 ±0.01 

0.86 

7.11 

x=0.1 OD 

0.657 ±0.01 

0.81 

7.06 

a;=0.1 MD 

0.561 ±0.02 

0.69 

6.94 

1=0.1 UD 

0.055 ±0.02 

0.07 

6.32 

1=0.4 UD 

0.067 ±0.01 

0.09 

6.34 


the OD samples, combined with the phase diagram for 
powders (see Fig. we obtain j/=7.06 and 7.II for 
the a:=0.4 and x=0.1 samples respectively. This is near 
the top, but slightly to the left of the peak of the su¬ 
perconducting domes. We then estimate the amount of 
holes using the relation h = y — 6.25^. Using that as a 
reference, and assuming the B/{A + B) area ratios are 
proportional to h, we can estimate h and y of the UD 
and MD samples. A summary of the intensity ratios, 
estimated h, and estimated y for the various samples is 
tabulated in tablelH We note that yuD — 6.32-6.34, plac¬ 
ing it well into the antiferromagnetic long-range ordered 
phase (Fig. [T]). Likewise, ?/mz> — 6.94, which is consis¬ 


tent with the iodometric titration result of 6.92 for this 
sample. 

To further compare the relative hole densities, the nor¬ 
malized oxygen iF-edge spectra is plotted in Fig. [T^l It 
was measured for the a:=0.1 OD, a::=0.4 OD, and x=0.1 
MD samples. The effect of the holes may be seen by in¬ 
spection of the positions of the low-energy peak of the O 
AT-edge spectra. Shifts in this oxygen “pre-edge” energy 
track the shift in Fermi level with hole dopin g^ This 
shift is a direct consequence of the filling (or emptying) 
of the bands. From Fig. [T^l the low-energy oxygen K- 
edges overlap almost exactly for the x = 0.4 and x = 0.1 
OD samples. In contrast, the edge of the x = 0.1 MD 
spectrum shifts by about 0.07 eV. Based on the result 
shown for YBCO in Ref. 52, the shift would correspond 
to a change in doping of 5y ~ 0.20. This is of the same 
order of magnitude as 5y ~ 0.12 between the OD and 
MD samples in table HI The almost overlapping edges for 
the a;=0.1 and a;=0.4 OD samples is therefore a second 
confirmation of identical number of holes, and further¬ 
more indicates that the amount holes in the plane layer 
are the same. 


(a) 

\ x=0.1 OD 







FIG. 13: (Golor Online) Fitting of the background-subtracted 
XAS spectra to three Lorentzians for the optimally doped 
samples for (a) 2;=0.1 and (b) a;=0.4. 
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